Broadband transmission properties of multilayered structures 
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The formalism of the scattering matrix is applied to describe the transmission properties of multilayered 
structures with deep variations of the refractive index and arbitrary arrangements of the layers. We show that 
there is an exact analytical formula for the transmission spectrum, which is valid for the full spectral range 
and which contains only a limited number of parameters for structures satisfying the quarter-wave condition. 
These parameters are related to the poles of the scattering matrix, and we present an efficient algorithm to 
find them, which is based on considering the ray propagation inside the structure and subsequent application 
of the harmonic inversion technique. These results are significant to analyze the reshaping of ultrashort pulses 
in multilayered structures. 

OCIS codes: 290.5825, 230.4170, 310.6860, 320.5540. 



Multilayered structures with periodic arrangement of the 
layers represent the simplest example of photonic crys- 
tals. Such structures are widely used as distributed re- 
flectors, spectral filters and can be applied for compres- 
sion or reshaping of ultrashort pulses [1]. A lot of at- 
tention has been paid to the question of how to mod- 
ify the properties of multilayered structures by intro- 
ducing artificial defects. As a more general case, struc- 
tures with quasiperiodic and deterministically aperiodic 
arrangements of the layers were considered [2] . The stan- 
dard technique to compute the transmission spectrum of 
multilayered structures is based on the transfer matrix 
method. However, being a strictly numerical method, it 
does not provide a proper understanding of the trans- 
mission properties from the physical point of view. 

In this paper, we show that there is an exact analytical 
formula for the transmission spectrum which is applica- 
ble for multilayered structures with a deep variation of 
the refractive index and arbitrary arrangement of the 
layers. It is valid for the full spectral range and contains 
only a limited number of parameters for structures satis- 
fying the quarter-wave condition. The knowledge of the 
precise analytical formula for the broad band transmis- 
sion allows one not only to compute easily many impor- 
tant characteristics such as the group velocity dispersion 
or the photonic density of modes, but also to analyze the 
propagation of arbitrary signals through the structure in 
the time domain. 

Any multilayered structure can be considered as a 
black box, the input and output from which are related 
by the scattering matrix. It is known that the poles and 
zeros of the scattering matrix determine the properties 
of a system uniquely [3]. Taking into account that the 
components of the scattering matrix can be interpreted 
as refiection and transmission coefficients, one can de- 
scribe the transmission spectrum T{ljo) of an arbitrary 
multilayered structure by the following formula 
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where ap is the strength of the resonance at cjp, and 
the sum is taken over all resonances. The resonances 
are associated with poles in the transmission spectrum 
T~^{ujp) = 0. For media without gain all poles are lo- 
cated in the lower part of the complex frequency plane 
lm{ujp) < 0. The strength of the resonances can be found 
as ap = [dT-\uj)/duj\-i^^. 

In general, the number of resonances is infinite, but 
due to the time-reversal symmetry T(— cj*) = T*(co'), 
which means that the poles exist in pairs {cjp, — cj*} and 
that their strengths are related as a{—cu*) = —a*{cjp). 
Those multilayered structures, which satisfy the quarter- 
wave condition rimdm = Aqw/4 at the wavelength Aqw for 
each layer m with the refractive index rim and thickness 
dm^ have the additional symmetry T(co'+4/cco'qw) = T{uj)^ 
where k is an arbitrary integer, cjqw is the frequency 
at which the quarter- wave condition holds {cu^y^/c = 
27r/Aqw). Therefore, there is only a limited number of 
independent poles P in the interval (— 2co'qw, 2co'qw], and 
the positions of all the others can be found by adding 
4/ca;qw. 

The inverse Fourier transform of Eq. (1) gives the re- 
sponse function G{t) of the structure to the excitation 
in the form of the Dirac delta function S{t) 
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where h{t) is the Heaviside step function, and the 
periodicity of the poles for the quarter-wave struc- 
tures was used explicitly. Applying the following 
property of the Dirac comb ^^=_^^{t — nT) = 
J^^_oo exp[— iA:(27r/r)t] leads to 

oo 

G{t) = J2anS{t-nT), (3) 

where r = 7r/(2u;qw) = Umdm/c has the meaning of 
time necessary for a signal to go through a quarter- 
wave layer. The coeflicients in Eq. (3) are defined as 
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Fig. 1. (Color online) (a, d) Transmission spectrum of a quarter- wave structure described by a binary sequence. The letters 
"A" and "B" correspond to layers with the refractive indices 1.55 and 2.3, respectively, (b, e) Positions of the poles in the 
complex frequency plane. In contrast to (a, d), where lm{(jj) = and |T(c(;)| < 1, transmission goes to infinity in the vicinity 
of poles, and the gray areas display regions where it exceeds a few fixed values, (c, f) Group delay as a function of frequency. 



an = —"^p^iiapTexp^—inujpr), and thus the response 
function can be viewed as a sum of exponentially decay- 
ing modes, which are sampled over the discrete intervals 
r. It is worth noting that all coefficients are real be- 
cause G{t) is real by definition. Moreover, the signal at 
the output cannot appear immediately, and = for 
n < M, where M is the total number of layers in the 
structure. These causality relations emphasize that <jp 
are not independent parameters and are related to uOp. 
The Fourier transform of Eq. (3) gives 
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which shows that are equal to the Fourier compo- 
nents of the transmission spectrum. Taking the sum of 
the geometric progression in Eq. (4) leads to 
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Therefore, the transmission through the quarter- wave 
structures can be described by a formula which does not 
involve infinite series and has only a limited number of 
parameters. We checked that it gives exactly the same re- 
sults as the transfer matrix method over the full spectral 
range. Two examples are shown in Fig. 1, which corre- 
spond to multilayered structures based on the Fibonacci 
sequence of the 7th order (on the left) and a periodic 
sequence with a defect (on the right). 

The periodicity of the poles could be used directly to 
obtain a few other formulas for the transmission. The 
repeated poles of equal strength can be constructed by 
using the expansion of cotangent into partial fractions [4] 
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Alternatively, one can operate with multiple zeros rather 



than poles and to rewrite the transmission as a product 
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The formulas (5)-(7), are fully equivalent to each other. 
For example, to derive Eq. (6) from Eq. (5), it is sufficient 
to notice that due to causality ao = cr^ = 0. 

To analyze the propagation of ultrashort pulses 
through the structure, it is important to compute the 
group delay as a function of frequency [5] . By separating 
the amplitude and phase in the transmission spectrum 
T{uj) = \T{uj) \ ex.-p{i(p), the group delay can be defined as 
rD(co') = dcp/duj = lm[T' {uj)/T{uj)]. An explicit formula 
for it can be obtained particularly easy from Eq. (7) 
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Since the integral J_^^^cot{z — Zp)dz = —iir for any 
Zp in the lower part of the complex plane Im[zp] < 0, 
the averaged variations of the group delay do not de- 
pend on the particular arrangement of the layers in the 
structure J^2^ rD(a;)da; = ttP. On the other hand, for 
a homogeneous slab consisting of M quarter-wave lay- 
ers {dip/duj)duj = 27rM. Therefore, P = 2M, or 
the number of resonances in the interval (— 2a;qw, 2a;qw] 
equals twice the number of layers [cf. Fig. 1]. 

The concept of group delay is closely related to other 
characteristics of multilayered structures such as the 
group velocity or the density of modes. This means 
that similar analytical formulas can be derived for them. 
For instance, the group velocity (or traversal velocity) 
can be defined as Vgj^ = L/tj:,^ where L is the to- 
tal length of the structure. It can be also rewritten as 
/c = Mr /{rB{n))^ where an averaged refractive index 
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)/L was introduced. This makes eas- 
ier the comparison of the data for the group delay shown 
in Fig. 1 with the group velocity. 
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Fig. 2. (Color online) Gaussian pulses before (solid) and 
after (dashed) propagation through a multilayered structure. 



In general, the transmitted signal /out(^) can be found 
as a convolution of G{t) with the input signal fin{t): 



fout{t)= G{t - t')fin{t')dt' = Y,anfin{t 
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This series converges very rapidly because the coeffi- 
cients an decay exponentially with increasing n. The 
simplest reshaping of the signals takes place when only 
one resonance is excited. As an example, we consider 
a well isolated resonance at ujp ~ ujq^ that exists for 
the structure shown in Fig. 1(d). It is worth noting 
that the group delay experienced by the pulses with 
the central frequency close to Re[cjp] can be estimated 
as rD(Re[co'p])/r ~ — 2/(7rIm[cjp/cjqw]), and it coincides 
with the decay time of the resonance in the absence of 
excitation. Which effect will dominate depends on the 
duration of the incident pulse [see Fig. 2]. The exponen- 
tial stretching is more pronounced for shorter pulses, and 
it can be accompanied with beating if several resonances 
are located close to each other [6]. 

Although the positions of the poles can be found by it- 
erations using the transfer matrix method, we developed 
a more reliable and efficient algorithm which is based on 
considering the ray propagation inside the structure. 

Initially, the amplitudes of forward and backward 
moving rays in all layers of the structure and on the 
boundaries are set to zero = OforO<m<M + l 
[see Fig. 3]. Then, a short probe signal is launched from 
the left, which corresponds to setting = 1 only for 
the first moment of time t = 0. Afterwards, a time 
marching scheme is applied to collect the outgoing sig- 
nals on the right boundary at discrete moments of time 
t = nr. For odd (even) time steps n, one should apply 
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Fig. 3. (Color online) Propagation of rays inside a quarter- 
wave multilayered structure. 



[«m^ ^m+i]^ Y^,m+i[Q^m^ ^m+iV with iudcx Tu run- 
ning over even (odd) numbers. The matrix Y^^^^ consists 
of the Fresnel coefficients which describe the scattering 
of rays at the interface between media with the refractive 
indices Uu and riy 
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Since the output signal on the right can be repre- 
sented as a sum of exponentially decaying modes [AJ^ = 
— ^^^-L icTpT exp(— mcjpr), the unknown amplitudes (jp 
and resonant frequencies uOp can be found as a solution 
of the following eigenvalue problem [7] 



HiVp = exp(-ia;pr)HoVp 



(11) 



where the matrix elements of auxiliary Hamiltonians H/g 
are defined as [H/c]^^^ = [A]^+^+/e. This harmonic in- 
version technique is very efficient because it requires the 
knowledge of only 2P elements in the vector A. The poles 
uOp are directly related to the eigenvalues of Eq. (11), and 
their strength <jp can be determined after a proper nor- 
malization of the eigenvectors Vp 
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In conclusion, we showed that the transmission spec- 
trum of multilayered structures can be described by ana- 
lytical formulas regardless of the specific arrangement of 
the layers. These formulas use the position of resonances 
as parameters, and we presented an efficient algorithm 
to find them. This establishes a new approach to ana- 
lyze the transmission properties of multilayers. It can be 
useful not only for an advanced reshaping of ultrashort 
pulses but also for the improvement of nonlinear effects 
such as self-pulsing or nonreciprocal transmission [8,9]. 
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